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ABSTRACT 

6 ■ 

■ We present the results of detailed spectral studies of the ultra-compact low mass 

X-ray binary (LMXB) 4U 1820-30 carried out with the Rossi X-ray Timing Explorer 

> ! (RXTE) during 1996-7. 4U 1820-30 is an "atoll" source X-ray burster (XRB) located 

in the globular cluster NGC 6624. It is known to have an 11 minute binary period 
and a ~ 176 day modulation in its 2-12 keV flux. Observations were made with the 
PCA and HEXTE instruments on RXTE at roughly one-month intervals to sample this 
long-term period and study flux-related spectral changes. We find that the 176-day 
period corresponds to normal atoll source motion in the color-color diagram, indicating 
that the long period is due to mass accretion rate changes. There are clear correlations 
between our fitted spectral parameters and both the broad-band (2-50 keV) flux and 
the position in the color-color diagram, as described by the parameter S a introduced by 
Mendez et al. (1999). The hard X-ray tail becomes more prominent at low flux levels as 
the source moves into the island state. In addition, we find a strong correlation between 
the position in the color-color diagram and the frequencies of the kilohertz quasi-periodic 
oscillations (kHz QPOs) reported by Zhang et al. (1998). This lends further support 
to the notion that evidence for the last stable orbit in the accretion disk of 4U 1820- 
30 has been observed. For a model consisting of Comptonization of cool photons by 
hot electrons plus an additional blackbody component, we report an abrupt change in 
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the spectral parameters at the same accretion rate at which the kHz QPOs disappear, 
suggesting a change in the accretion flow that affects both the spectrum and the timing 
properties. For a model consisting of a multicolor disk blackbody plus a cut-off power 
law, we find that the inner disk radius reaches a minimum at the same accretion rate 
at which the kHz QPO frequency saturates, as expected if the disk reaches the last 
stable orbit. Both models face theoretical and observational problems when interpreted 
physically for this system. 

Subject headings: accretion, accretion disks — stars: individual (4U 1820-30) — stars:neutron 
— X-rays:stars 

1. Introduction 

The broad-band spectral sensitivities of two X-ray astronomy spacecraft, the Rossi X-ray 
Timing Explorer (RXTE, Bradt, Rothschild, & Swank 1993) and BeppoSAX (Boella et al. 1997), 
combined with the high timing resolution of RXTE, have allowed major advances in recent years in 
the study of Low Mass X-ray Binaries (LMXBs). Progress has been particularly rapid in the field of 
X-ray bursters (XRBs), binary systems containing a low-magnetic field neutron star (NS) accreting 
matter from a late-type companion star via Roche lobe overflow. It had become apparent by the 
mid 1990's that NS systems were capable of producing hard X-ray (> 20 keV) emission, as indicated 
by BATSE and SIGMA detections of several XRB at ~ 100 keV (Barret & Vedrenne 1994, Tavani 
& Barret 1997 and references therein). This eliminated hard X-ray radiation as a unique identifier 
of black hole (BH) systems, but the limited sensitivity and spectral resolution of these instruments, 
together with the lack of simultaneous soft X-ray observations, did not allow the spectral shape to be 
studied in detail or possible differences between NS and BH hard X-ray properties to be determined. 
The two new satellites are revolutionary due to their broad spectral coverage: the Proportional 
Counter Array (PCA) and High Energy X-ray Timing Experiment (HEXTE) instruments on RXTE 
cover the range 2-150 keV, while the suite of narrow-field instruments (NFIs) on BeppoSAX combine 
to cover 0.1-200 keV. Many XRB have now been observed with these experiments, allowing detailed 
modeling of the complete X-ray to hard X-ray spectra with a variety of complex models (e.g., Church 
et al. 1998; Guainazzi et al. 1998; in 't Zand et al. 1999; Olive et al. 1999, Piraino et al. 1999a; 
Piraino et al. 1999b; Barret et al. 1999, 2000). 

In addition to broad spectral coverage, RXTE offers unprecedented timing resolution (~ l//s) 
due to the large collecting area of the PCA. This has led directly to the discovery of quasi-periodic 
oscillations (QPOs) with frequencies near 1000 Hz (kHz QPOs) in the X-ray timing power spectra 
of many NS systems (Strohmayer et al. 1996; van der Klis et al. 1996). Since kilohertz frequencies 
are naturally associated with orbital timescales near a NS, it is generally accepted that the kHz 
QPOs are produced in the inner accretion disk and that the frequency is related to the inner 
disk radius. If this is true, the kHz QPOs are a potentially valuable probe of strong field gravity 
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(Kaaret & Ford 1997). In particular, the relationship between QPO frequency and accretion rate 
should provide a means to test for the existence of the marginally stable orbit predicted by general 
relativity (Kaaret, Ford, & Chen 1997; Zhang, Strohmayer, & Swank 1997). 

In this paper we present the results of RXTE observations of 4U 1820-30, an XRB with a 
particularly colorful history. 4U 1820-30 is located in the core of the globular cluster NGC 6624 
and was one of the first globular cluster X-ray sources discovered (Giacconi et al. 1974). The 
distance to NGC 6624, 6.4 ± 0.6 kpc, was derived by Vacca et al. (1986) from a study of horizontal 
branch stars. Thus the distance to 4U 1820-30 is among the better known of all X-ray sources. 
It is the brightest globular cluster source, and was the first identified X-ray burster (Grindlay et 
al. 1976). The source has been observed to shift between a high and a low state with a period of 
about 176 days (Priedhorsky & Terrell 1984); in the All Sky Monitor (ASM) on RXTE the 2-12 
keV flux varies by typically a factor of ~ 3 and as much as a factor of ~ 6 between the two states 
(see Figure 1). Also observed is an 11 minute (685 s) period believed to be the orbital period, the 
shortest known orbital period of all cosmic binary systems (Stella, Priedhorsky, & White 1987). 
Such a short period suggests that the mass-donating star is a helium dwarf with a mass of 0.06-0.08 
Mq (Rappaport et al. 1987). Two simultaneous kHz QPOs were detected in the persistent flux 
from 4U 1820-30 by Smale, Zhang, & White (1997). 

The goal of the present investigation is to observe the broad-band (2-150 keV) spectrum of 4U 
1820-30 with the PCA and HEXTE instruments on RXTE and to study the relationships between 
spectral parameters, kHz QPO frequencies, and accretion rate. The steady 176-day variation of the 
flux from 4U 1820-30 is most likely due to modulation of the accretion flow, since X-ray bursts are 
only observed in the low state (Stella, Priedhorsky, & White 1987), and so this source provides a 
good laboratory for the study of accretion-related spectral changes. We have observed 4U 1820-30 
over the course of a year at roughly monthly intervals to sample the entire 176-day period. This 
same data set was used to study timing properties by Zhang et al. (1998). In Section 2 we describe 
previous spectral and timing observations of 4U 1820-30. In Section 3 we describe the observations 
and our analysis procedure. The results are presented in Section 4, discussed in Section 5, and we 
state our conclusions in Section 6. 

2. Previous Observations of 4U 1820-30 

As evidenced by its long history, 4U 1820-30 is an extremely well-studied system. In this 
section we briefly review previous spectral and timing observations of the source. 

2.1. Spectral Observations 

4U 1820-30 has previously been observed by ANS, EXOSAT, Ginga, ASCA, Einstein, and Bep- 
poSAX. A wide variety of models have been used to fit the recorded spectra, but the most successful 
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have been two-component models consisting of either optically thin thermal bremsstrahlung plus a 
blackbody (TB + BB), a power law with exponential cutoff (an analytical approximation of unsat- 
urated Comptonization) plus a blackbody (CPL + BB), or a more detailed Comptonization model 
developed by Sunyaev & Titarchuk (1980) plus a blackbody (CompST + BB, after the CompST 
model in XSPEC). Models based on optically thin thermal bremsstrahlung are widely believed to 
be unphysical for XRB, since the emission measures typically derived (~ 10 60 cm" 3 ) imply, for 
small optical depths, that the plasma cloud has a radius ~ 10 10 cm, too large to be heated to 
several keV by the neutron star (White et al. 1986; Christian & Swank 1997). We therefore restrict 
our attention to the CPL and CompST (with or without a BB) model results. The important 
parameters for the CPL + BB model are the blackbody temperature kTsB, the power law photon 
index a, and the exponential cutoff energy E c . For the CompST + BB model, the parameters of 
interest are ¥Tbb and the temperature kT e and optical depth r of the Comptonizing electron cloud, 
assumed here to be spherical. Photoelectric absorption at low energies by the interstellar medium 
must be included as well, given by the equivalent hydrogen column density Njj (found using the 
absorption cross sections of Morrison &i McCammon(1983)). Finally, a gaussian line feature at 
~ 6.7 keV is often required, representing a blend of Ka emission lines from helium-like ions of 
iron (Fe xxv or XXVI; White et al. 1986; Hirano et al. 1987). A summary of previously-reported 
spectral fits using these models is given in Table 1 and shown graphically in Figure 10; we discuss 
them briefly here. 

Intensity-related spectral changes were first reported for 4U 1820-30 by Parsignault & Grindlay 
(1978) using ANS/EXX data and a power law fit. They found that the photon index a varied 
between ~ 2 at high count rates and ~ 1.4 at low count rates, indicative of spectral hardening at 
low luminosity. More complex models could be used after the EXOSAT MK instrument observed 
4U 1820-30 between 1984 and 1985 (Stella, White, & Priedhorsky 1987), finding a high state 
luminosity (1-30 keV) of 6.0 x 10 37 ergs s _1 and a low state luminosity of 2.0 x 10 37 ergs s _1 . Fits 
with the CPL + BB model (Stella, White, & Priedhorsky 1987) gave low state values of a = 2.5, 
E c > 30 keV, and UTbb = 2.3 keV, while in the high state a decreased to 1.7, E c decreased to ~ 12 
keV, and UTbb fell to 2.0 keV. Note that the inclusion of an exponential cutoff has changed the 
intensity-related behavior of a compared to the ANS results; the softening of the spectrum is now 
described by the lower cutoff energy while a describes the spectrum below this. White, Stella, & 
Parmar (1988) fitted the CompST + BB model to the EXOSAT data and found that between the 
low and high states VTbb decreased from 2.06 to 1.61 keV, kT e fell slightly from 3.5 to 3.3 keV, 
and r decreased from 13.2 to 11.9. For both models the ratio of blackbody to "hard component" 
flux Lbb/Lh decreased from 50% to 25% as the luminosity increased. In a similar energy band, 
the Ginga LACs observed 4U 1820-30 at two different intensity levels in 1987 May (Ercan et al. 
1993), but the power law was steeper with a lower cutoff than in the EXOSAT spectrum (Table 1). 
The ratio Lbb/Lh was similar however. 

Lower energy coverage was provided by ASCA/GIS and Einstein observations. ASCA observed 
4U 1820-30 in 1993 in the low state and the spectrum (0.6-11 keV) was fitted with the CompST + 
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BB model (Smale et al. 1994). CompST parameters were similar to those found by EXOSAT, but 
the blackbody temperature was substantially lower {kT BB = 0.76 ± 0.02 keV). Christian & Swank 
(1997) report on an Einstein (SSS + MPC) observation (0.5-20 keV) of 4U 1820-30 from 1978. 
The luminosity was 5.5 x 10 37 ergs s _1 , indicative of the high state, and yet fitting the CompST 
+ BB model gave paramters very similar to those from ASCA (Table 1). The blackbody flux was 
only 6% of the CompST flux. 

The first observations of 4U 1820-30 extending above ~ 30 keV were made by the NFIs (0.1- 
200 keV) on BeppoSAX in 1998 (Piraino et al. 1999a; Kaaret et al. 1999). These were the first 
truly broad-band observations, though they covered only a narrow range in luminosity (Table 1). 
A CPL + BB fit to the 0.3-40 keV spectrum gave kT BB = 0.47 keV, a = 0.55, and E c = 4.5 keV, 
while the CompST + BB model yielded kT BB = 0.46-0.66 keV, kT e = 2.83 keV, and r = 13.7. 
Thus the instruments with response below 1 keV consistently find low values of kT BB . The entire 
range of luminosities observed so far with all instruments has yielded very similar values for the 
parameters of the CompST + BB model. 

The BeppoSAX Phoswich Detection System (PDS) failed to detect 4U 1820-30 above 40 keV, 
consistent with the failure of BATSE to detect the source (20-100 keV) over the first four years of 
the CGRO mission (Bloser et al. 1996). Thus 4U 1820-30 is not a member of the class of XRB 
detected at ~ 100 keV. 

2.2. Kilohertz QPOs 

Two simultaneous kHz QPOs were discovered in the persistent emission of 4U 1820-30 by 
Smale, Zhang, & White (1997) using RXTE/PCA data from 1996 October. The characteristics 
of the QPOs were similar to those found in other XRB: their centroid frequencies were correlated 
with the 2-60 keV PCA count rate, but the frequency difference remained constant at 275 ± 8 Hz. 
The fractional rms amplitudes increased dramatically with energy. Zhang et al. (1998) studied 
monthly PCA observations of 4U 1820-30 (the same data set reported on here; see Section 3) and 
found that the frequencies of the QPOs are only correlated with PCA count rate below a critical 
value. Above this, the QPO frequencies remain constant while the count rate increases. Saturation 
of QPO frequency at high mass accretion rates is an expected signature of the marginally stable 
orbit (Miller, Lamb, & Psaltis 1998; Kaaret, Ford, & Chen 1997); however, Mendez et al. (1999) 
have shown that PCA count rate is not a good indicator of accretion rate. Kaaret et al. (1999) 
used more reliable estimators of mass accretion rate, broad-band flux and hard color, to show that 
the frequency saturation is robust and thus strong evidence that the inner disk in 4U 1820-30 does 
reach the marginally stable orbit at high accretion rates. 
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3. Observations and Analysis 

3.1. RXTE Observations and Data Reduction 

The PCA instrument on RXTE is made up of five Xenon proportional counter units (PCUs, 
numbered 0-4) sensitive from 2-60 keV with a total area of about 6500 cm 2 (Bradt, Rothschild, & 
Swank 1993). The HEXTE instrument consists of two clusters of four NaI(Tl)/CsI(Na) phoswich 
scintillation detectors sensitive from 15-200 keV with an effective area of 1600 cm 2 ; these clusters 
rock on and off source to obtain background measurements (Rothschild et al. 1998). The ASM is 
comprised of three Scanning Shadow Cameras (SSCs) with one-dimensional slit masks and 6° x 90° 
fields of view that monitor the fluxes of bright X-ray sources several times a day in the 2-12 keV 
band (Levine et al. 1996). 

Sixteen joint PCA/HEXTE pointed observations of 4U 1820-30 were performed at roughly- 
monthly intervals between 1997 February 9 and 1997 September 10. In addition, seven observations 
made between 1996 October 26 and 31 were included (Smale et al. 1997). The RXTE/ ASM 2- 
12 keV light curve of 4U 1820-30 during this period is shown in Figure 1, with the times of our 
pointed observations marked. Each point represents the one-day average of the ASM count rate. 
The 176-day modulation is clearly visible and is well-sampled by the pointed observations. The 
ASM count rate varies from ~ 5 to ~ 30 counts s _1 during this time, or ~ 67 mCrab to ~ 400 
mCrab. Assuming a Crab- like spectrum and Njj = 0.3 x 10 22 cm~ 2 (see below), these low and 
high state values correspond to 1-20 keV luminosities of ~ 1.2 x 10 37 ergs s _1 and ~ 7.0 x 10 37 
ergs s _1 , respectively, consistent with the range of luminosities observed previously (Table 1). The 
dates, starting times in seconds of Mission Elapsed Time (MET = seconds since 1994 Jan 1, OhOmOs 
UTC), and mean PCA count rates (2-16 keV) for all 23 observations are listed in Table 2. 

The data were divided into the same 90 segments (persistent emission only, typical length 
~ 3000 seconds) used by Zhang et al. (1998) and reduced using the standard RXTE analysis tools 
contained in FTOOLS 4.2. For our color-color diagram and PCA spectral analysis we used the 
"Standard 2" data, which provide 128 energy channels between 2 and 100 keV with 16-second 
time resolution. Spectra were extracted with saextrct v4.0b. Backgrounds were estimated with 
pcabackest v2.1b (October 1998), and the v2.2.1 response matrices (January 1998) were used. For 
the HEXTE spectral analysis we used event mode (E_8us_256_DXlF) data with the 1997 March 
response matrices. Event mode data were used instead of the standard archive mode data because 
the 1996 archive data files were corrupted; this problem has been corrected in the reprocessed 
distribution of 1996 data. Spectra were extracted with saextrct and rebinned by a factor of two 
using grppha to provide the same energy bins as standard archive data below 60 keV. Deadtime 
corrections were computed using hxtdead vO.0.1. 
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3.2. Color-Color Diagram and Parameterization of Accretion Rate 

A color-color diagram (CCD) was produced from the PCA data with one point for each of the 
90 data segments. Only PCUs 0, 1, and 2 were on continuously during all observations, and so 
only data from these PCUs are included here. The soft and hard colors are defined as the ratios of 
background-subtracted PCA count rates in the bands 3.5-6.4 keV and 2.0-3.5 keV, and 9.7-16.0 
keV and 6.4-9.7 keV, respectively. The CCD is shown in Figure 2. We observe 4U 1820-30 in all 
three of the states associated with "atoll" sources (Hasinger &: van der Klis 1989): the island state, 
lower banana branch, and upper banana branch. This confirms that the 176-day modulation is in 
fact due to variation in the mass accretion rate. Figure 2 shows that we have adequately sampled 
the long period and observed the source in all of its possible states. Closed symbols indicate data 
segments in which kHz QPOs were discovered by Zhang et al. (1998), while filled symbols indicate 
data without kHz QPOs. 

It is generally believed that in atoll sources the mass accretion rate M increases monotonically 
as the source moves from the island state, through the lower banana, and into the upper banana 
branch (Hasinger & van der Klis 1989). Mendez et al. (1999) defined the parameter S a to measure 
the position of an atoll source within its CCD and thus parameterize its accretion rate. We follow 
Mendez et al. (1999) and approximate the shape of the CCD track with a spline, as shown in 
Figure 2. S a is then defined as the distance along this curve, as measured from the island state. 
We arbitrarily set S a = 1.0 at (2.378,0.650) and S a = 2.0 at (2.671,0.597) in the CCD. Each data 
point is assigned the value of S a of the point closest to it on the curve. The error in S a is estimated 
by simply projecting the error bars of the data points onto the curve. Thus we may use in 
addition to the source luminosity to study accretion rate-related changes in spectral parameters or 
kHz QPO behavior. The advantage in using S a instead of the hard color (Kaaret et al. 1999) to 
track the dependence of QPO frequency on accretion rate is that S a is more sensitive to changes 
of state along the nearly horizontal banana branch, which is where the break in QPO frequency 
observed by Zhang et al. (1998) occurs. 



3.3. Spectral Fitting 

The spectral analysis was performed using XSPEC vlO.O (Arnaud 1996). Spectra for each 
PCU and HEXTE cluster were reduced separately and combined within XSPEC. Due to systematic 
uncertainties in the response matrices, only HEXTE data above 20 keV was used. Several authors 
have written on systematic effects in the PCA and their effect on spectral fitting (e.g. Sobczak 
et al. 1999; Barret et al. 2000; Wilms et al. 1999; Tomsick et al. 1999); we performed similar 
investigations and reached similar conclusions. We extracted spectra for each PCU from archival 
Crab observations spanning the year of our observations (1996 Oct 27, 1997 Mar 22, 1997 Jul 26, 
and 1997 Dec 15) and fit them with an absorbed power law. Initially we used response matrices 
generated by pcarmfv3.5, but we got very poor fits with large systematic deviations in the residuals, 
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especially around 5 keV. The residuals were far better using the v2.2.1 response matrices for all 
four Crab observations, and so we used these for all our 4U 1820-30 fits. The residuals still showed 
large systematic effects below 2.5 keV and above 20 keV, and so we restricted our fits to this energy 
range. We also found, as have other authors, that the residuals in PCUs 2 and 3 are considerably 
worse than those in PCUs 0, 1, and 4. For this reason we use only PCUs 0, 1, and 4 in our 
fits. PCU 4 was on during all observations except 1997 May 28. The relative normalizations of 
the three PCUs and two HEXTE clusters were left as free parameters, since there still exist large 
uncertainties in the relative flux normalizations. In our Crab fits, all three PCUs find power law 
normalizations 20%-30% higher than the accepted value of ~ 10 photons cm~ 2 s _1 keV^ 1 at 1 
keV (e.g. Jung 1989). Tomsick et al. (1999) have reported similar results. All spectral model 
normalizations reported here are obtained from PCU 0; in Tables 3- 5 (and also Figure 5) we 
report the flux values found directly by PCU 0, while when comparing luminosities with previous 
measurements (Figure 10) we reduce the PCA values by 15% (see Section 5). A systematic error 
of 1% was added to all PCA channels using grppha. In Figure 3 we show two representative PCA 
and HEXTE spectra of 4U 1820-30, one from the low island state and one from the high banana 
state. The hardening of the spectrum at low luminosity is clear. In the island state the source is 
detected by HEXTE only up to 50 keV, confirming the lack of a detection above this by BeppoSAX 
(Piraino et al. 1999) and BATSE (Bloser et al. 1996). 

For our spectral analysis data segments were grouped according to their position in the CCD 
as follows: Spectra were extracted for each of the 90 data segments. The four island state segments 
were combined into one spectrum, and the five uppermost banana state segments were each analyzed 
as separate spectra. In between, data segments were grouped into 26 bins of width 0.025 in S a and 
combined together, with an average value of S a and the total exposure time computed for each bin. 
Thus we had a total of 32 spectra to fit (see Tables 3- 5). 

All spectra were initially fit with two models: the CPL model described in Section 2, and 
the recently-developed Comptonization model of Titarchuk (1994), called the CompTT model in 
XSPEC. This model improves on the CompST model by including relativistic effects and extending 
the theory to work in both the optically thin and thick cases. In addition to the Comptonizing 
electron temperature kT e and optical depth r, an additional parameter is the temperature of the 
cool seed photons kT\y, assumed to follow a Wien law. We consider only a spherical geometry for 
comparison with previous results. In nearly all cases the model fits were improved by the addition 
of a blackbody component. We also included a gaussian line component at ~ 6.7 keV in all spectra 
and a smeared absorption edge in the island state, as described in more detail in Section 4.2. 

Some XRB have recently been fit with a multi-color disk blackbody (DBB; Mitsuda et al. 
1984) instead of a simple single-temperature BB to account for the soft excess (Guainazzi et al. 
1999; Barret et al. 1999). We also fit all our spectra with a model consisting of a CPL plus this 
disk blackbody (CPL + DBB model). The DBB model parameters are the color temperature of 
the inner disk /cTj n and the projected inner disk radius Ri n V cos 9, where 9 is the inclination of the 
system. It is expected, however, that the true spectrum emitted by the accretion disk will be that of 
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a "diluted blackbody" due to the effects of electron scattering, and the DBB parameters must then 
be modified by a spectral hardening factor / (Ebisawa et al. 1994). The effective temperature and 
inner radius are kT e ff = kTi n / f and R e ff = f 2 Rin- Shimura & Takahara (1995) have found that 
/= 1.7 is a good approximation for accretion parameters appropriate for an XRB. We attempted 
also to fit a CompTT + DBB model, but the data were unable to constrain kT e and r for this 
combination. 

The insensitivity of the PCA below 2 keV did not allow us to determine the hydrogen column 
density Nh from our spectral fits. It was therefore necessary to freeze Nh at previously-determined 
values for all spectra. Optical reddening for NGC 6624 has been measured by Liller & Carney (1978) 
and Peterson (1993) to be Eb~v ~ 0.27, which, using the canonical relationship Ay = 3Eb-Vi 
yields an optical extinction Ay = 0.81. Predehl & Schmitt (1995) find that N H = 1.79 x 10 21 Ay 
cm~ 2 , and this gives a value of Nh ~ 1.4 x 10 21 cm -2 for NGC 6624. Previous X-ray observations 
of 4U 1820-30 with good low energy coverage (Einstein, BeppoSAX; Smale et al. (1994) do not 
report a value from ASCA) gave N H ~ 2.9 x 10 21 cm" 2 , and EXOSAT found N H ~ 4 x 10 21 cm" 2 
(see Table 1). Predehl & Schmitt (1995) list values of 2.1-3.6 xlO 21 cm" 2 for 4U 1820-30 from 
ROSAT data. These values are all in good agreement, and we therefore took a rough average and 
set Nh = 3.0 x 10 21 cm -2 and kept it frozen for all our spectral fits. 



4. Results 

4.1. Saturation of QPO Frequency 

The frequencies of the kHz QPOs reported by Zhang et al. (1998) are plotted as a function 
of S a , as derived here, in Figure 4. Open and filled plotting symbols indicate the lower and higher 
frequency QPOs, respectively, as defined by Zhang et al. (1998). There is a clear correlation of 
QPO frequency with S a below S a ~ 1.3, and then the frequency saturates even though S a continues 
to rise. This is an expected signature of the inner disk reaching the marginally stable orbit (Miller, 
Lamb, & Psaltis 1998). The QPOs disappear completely at S a ~ 1.45. These two "special" values 
of S a are indicated in Figure 2. If S a , the position of the source in the CCD, may indeed be taken 
to be an indicator of M , then Figure 4 provides strong evidence for the detection of the marginally 
stable orbit in 4U 1820-30. 

We have also investigated the possibility that the shape of the CCD varies with the 685-second 
orbital period. The CCD was reproduced with bins of 128 seconds, and each point was grouped 
into one of five phase bins. No significant change in the track described by the CCD was found, 
as the scatter in the points lay well within the error bars. In addition, we have confirmed that the 
shape of the frequency- S a plot in Figure 4 does not become smeared out when smaller time bins 
(and their associated larger error bars) are used. 
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4.2. Evolution of Spectral Parameters 

The 32 combined PCA/HEXTE spectra were first fit with the CPL + BB model with a gaussian 
line whose parameters were free to vary. Figure 5 shows the gaussian centroid energy E\ and width 
(a) as a function of the total luminosity in the 2-50 keV band spanned by both instruments (see 
also Table 3). In all but the low island state observation the gaussian parameters are consistent 
with constant values of E\ ~ 6.8 keV and a ~ 0.8 keV. These are reasonable values for an Fe xxv 
or xxvi Ka line reflected from the inner regions of a hot disk (White et al. 1986; Hirano et al. 
1987). The island state line energy is 6.53 keV, suggestive of reflection from less highly ionized gas, 
but the difference is only ~ 2a and it is difficult to see how the ionization could change so much 
for only a ~ 14% change in luminosity. We therefore froze E\ at 6.8 keV and a at 0.8 keV for our 
fits with the CompTT + BB model, since we found that the gaussian parameters were less well 
defined with this choice of continuum. The gaussian parameters were left free for the CPL + DBB 
fits, and no difference was found between the island and banana state line energies (Table 5). 

Since a broad iron line suggests reflection from a disk, we investigated the presence of other 
reflection components in our spectra. We first attempted to fit a smeared absorption edge (Ebisawa 
et al. 1994; smedge model in XSPEC) between 7 and 8 keV. Only in the island state spectrum could 
this component be fit with reasonable values. Although we would expect reflection to increase as 
the accretion rate and disk area increase, this model requires a power law continuum to fit against, 
and in all the rest of our observations the cutoff energy is too low, resulting in too "curvy" a 
spectral shape for reliable fitting of the edge with our data. The edge energies in the island state 
are 8.90+°;^ keV for the CPL + BB model and 7.73+J;^ keV for the CompTT + BB model, 
indicative again of ionized material. The smearing width was frozen at 10 keV. These edges are 
included in the island state fits shown in Tables 3 and 4. The smeared edge could not be fit with 
the CPL + DBB model. We attempted next to fit the complete XSPEC reflection model pexrav 
(Magdziarz & Zdziarski 1995) to the island state spectrum. This model is the combination of a 
cut-off power law with a Compton reflection component. The derived reflection scaling factor was 
consistent with zero for cos 9 = 0.45. We conclude that, other than the broad iron line, we have 
only marginal evidence for reflection from the disk, and only in the low state. 

The best-fit spectral parameters for the CPL + BB model are given in Table 3, the parameters 
for the CompTT + BB model are given in Table 4, and the parameters for the CPL + DBB model 
are given in Table 5. In all three Tables we give the average value of S a for the included data, and 
in Table 3 we list the total integration time for the PCA and HEXTE. The CPL + BB spectral 
parameters are as in Table 1, with the addition of the effective radius of the BB emitting surface 
Rbb and the equivalent width of the gaussian line. The luminosity and flux ratio Lbb/Lh are 
given for the 2-50 keV band covered by the PCA and HEXTE. (This is the raw luminosity derived 
directly from the fit.) All of the fits are acceptable, and in many cases the xt is so l° w as to indicate 
that the systematic errors have been overestimated. In Table 4 we give the best-fit temperature of 
the seed photons kTyy and the Comptonizing y-parameter, defined as y = 4kT e T 2 /m e c 2 . We also 
follow in 't Zand et al. (1999) and derive the effective Wien radius R\y of the seed photons: the 
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total flux in the CompTT component is corrected for energy gained by Comptonization by dividing 
by 1 + y, and this is set equal to the Wien flux from the surface of a sphere. The expression used 
is Rw = 3 x lt^D-^FcompTT/O- + y)/{kT\v) 2 km, where D is the distance in kpc, FcompTT is in 
erg cm" 2 s -1 , and kT\y is in keV. The derived values are listed in Table 4. The quality of the 
CompTT + BB fits is practically identical to that of the CPL + BB fits; based on the values of xl 
there is no reason to favor one model over the other. We note however that several of the CompTT 
+ BB fits have very large errors on the BB parameters. In Figure 6 we show two examples of 
spectra from PCUs 0, 1, and 4 and HEXTE clusters A and B fit with the CompTT + BB model. 
These are the same two spectra shown in Figure 3; now we show the raw counts with the folded 
model, the residuals, and the unfolded spectra with each component plotted separately. We note 
that the "wave" shape seen in the residuals of PCA spectra around 6 keV is still visible. In Table 5 
we include the fitted disk parameters kTi n and Ri n V cos 9 with the standard CPL and gaussian 
parameters. The fits are again all formally acceptable. We note, however, that in the island state 
the model is unable to reproduce the measured counts above 40 keV. 

In Figure 7 we show the spectral parameters of the CPL + BB and CompTT + BB models 
as a function of S a , and thus presumably of M . In both models it is clear that the parameters of 
the hard spectral component are correlated with the accretion rate, especially as the accretion rate 
become low, in the sense that the spectrum becomes harder as M decreases. The cutoff energy 
E c in the CPL model and kT e in the CompTT model both rise as the accretion rate drops. Note 
that when fitting a CPL model it is not a but E c that indicates the hardness of the spectrum, 
since a here is really an indication of the optical depth. This change is not gradual throughout 
the CCD, however; only as the source moves from the lower banana to the island state do the 
parameters change appreciably. The optical depth r of the CompTT model drops at lower M , 
and the power law index a (inversely related to r in the CPL approximation) rises. This general 
behavior is consistent with that seen in other sources (e.g., Barret et al. 2000; Christian & Swank 
1997). 

The most intriguing feature of Figure 7 is that the parameters of the CompTT + BB model 
show a distinct change at S a ~ 1.45. The ratio Lbb/Lh is inversely correlated with S a below this, 
but becomes constant above. The blackbody temperature UTbb drops from ~ 2.3 keV to ~ 1.6 
keV, and the blackbody radius Rbb begins to increase from ~ 1 km to ~ 2 km. The r versus S a 
curve also seems to roll over and become nearly constant at this value of S a . In fact, r and Lbb/Lh 
appear strongly anti-correlated below S a = 1.45. No such abrupt change is seen in the parameters 
Lbb/Lh, kTsB, Rbb, or a of the CPL + BB model. The cutoff parameters E c in the CPL model 
and kT e in the CompTT model do not seem to change at S a = 1.45; however, both appear to 
change from strong anti-correlations with S a to near-constant values at S a ~ 1.3. These "special" 
values of S a = 1.3 and S a = 1.45 are especially interesting because, as shown in Section 4.1 above, 
these are the same values at which the kHz QPOs saturate in frequency and disappear completely, 
respectively. Thus the accretion flow must be changing at the values of M corresponding to these 
values of S a in a way that affects both the energy spectrum and the frequencies of the kHz QPOs. 
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The fact that this correlation is seen only with the CompTT + BB model suggests that this model 
is the more accurate description of the spectrum. In Figure 8 we show the seed photon temperature 
kTyy and the Wien radius Ry/ described above as a function of S a . Although the error bars are 
large, there is an indication that the temperature and radius are anticorrelated below S a = 1.3 and 
constant above. This is precisely the behavior expected if the seed photons are coming from the 
inner part of the accretion disk, and the disk is contracting until it reaches the last stable orbit at 
S a = 1.3 as indicated by the kHz QPOs. This trend gives us confidence in the fitted values of kTyy, 
which must still be treated with some caution because the peak of the Wien spectrum at 3kT\y 
lies below the PCA energy range. The large values of the Wien radius also suggest that the seed 
photons originate in the disk. 

The parameters for the CPL + DBB model are shown in Figure 9. The cut-off power law 
components behave similarly to those in the CPL + BB model, though the values of a and E c 
are both considerably lower. This is because the DBB model prefers lower temperatures than the 
simple BB, and so the curvature in the spectrum at higher energies must now be described by a 
lower cutoff value. The behavior of the DBB component is of great interest: fcTj n rises and Ri n 
falls with increasing accretion rate until S a ~ 1.3, then both remain nearly constant (though there 
is some indication that kTi n continues to rise slowly). The fraction of the total luminosity in the 
DBB component increases with S a in the same manner. Thus it appears that the accretion disk 
inner radius grows smaller and hotter as the accretion rate increases until the last stable orbit 
is reached at S a = 1.3. This behavior of the DBB component is consistent with that expected 
from the behavior of the kHz QPOs. Thus this model is also compelling as an accurate physical 
description of the system. 



5. Discussion 

5.1. Comparison with Previous Spectral Measurements 

First we compare our spectral results for 4U 1820-30 to those previously reported using the 
CPL + BB and CompTT + BB models. In Figure 10 we show our RXTE spectral parameters 
together with the previously-reported parameters described in Section 2.1 and Table 1. Here we plot 
the spectral parameters as a function of the total 2-50 keV luminosity, assuming a distance of 6.4 
kpc. The luminosities of the previous measurements have been translated into the 2-50 keV band 
based on the reported parameters and fluxes. The integrated 2-50 keV flux in the archival Crab 
spectra (see Section 3.3) is ~ 15% higher than what is expected; therefore, the RXTE luminosities 
in Figure 10 have been divided by a factor of 1.15 for proper comparison to previous instruments. 
It is immediately obvious that there is considerable disagreement between different instruments. In 
general, the RXTE spectral parameters agree most closely with those found by EXOSAT and Ginga, 
and differ the most from those found by BeppoSAX. More specifically, the RXTE BB component 
parameters agree fairly well with those found by EXOSAT and Ginga for the CPL + BB model 
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and with those found by EX OS AT for the CompTT + BB model. In the CompTT + BB model, 
for example, EX OS AT also observed KTbb to decrease from > 2 keV to ~ 1.5 keV between the low 
and high states, while Rbb increased from ~ 1.7 km to ~ 3.7 km. The EXOSAT CPL parameters 
also agree well with our values in the low and high states, as do the Einstein values (Christian 
& Swank (1997) do not report BB values for the CPL model with Einstein). The BeppoSAX 
BB parameters and CPL spectral index a differ considerably from ours. All of the experiments 
that used a CompST model found Comptonization parameters that agree well with our CompTT 
parameters at high luminosity. In the low state, however, EXOSAT did not observe the increase 
in kT e and decrease in r that we see. It is strange that the EXOSAT low state observation could 
only place a lower limit of ~ 30 keV on E c for the CPL model, yet could still constrain such a 
low kT e ~ 3.5 keV for the CompST model. It is apparent that our island state observation is 
the first observation of 4U 1820-30 in the low state since EXOSAT, and it is puzzling that the 
Comptonization parameters differ so much. Further observations in the island state are needed, 
but the low duty cycle in this state seen in Figure 1 makes this difficult. 

It is apparent from Table 1 and Figure 10 that those experiments without sensitivity below 1 
keV, i.e. EXOSAT, Ginga, and RXTE, find consistently higher UTbb and lower Rbb than do those 
instruments with sensitivity below 1 keV, i.e. ASCA, Einstein, and BeppoSAX. The differences are 
important to the question of the origin of the soft component of the spectrum. The instruments 
with low energy sensitivity find Rbb <; 10 km, consistent with the radius of a neutron star. The 
other instruments, including RXTE, find Rbb ^ 4 km, implying that the blackbody component 
can come from only a limited region of the neutron star surface, such as an optically thick boundary 
layer. (An alternate explanation for apparently small blackbody radii for XRBs given by Rutledge 
et al. (1999), that a blackbody is not an appropriate model and fits with hydrogen atmosphere 
models give radii near 10 km, is not applicable to 4U 1820-30 since the accretion rate is too high.) 
We have attempted to fit our RXTE spectra with a BB component with kTs b < 1 keV, forcing the 
fit to find a new local minimum. For a CompTT + BB model of this type, VTbb remained nearly 
constant at ~ 0.6 keV while Rbb stayed at ~ 24 km, too large for a NS surface. The seed photon 
energy kT\y decreased from 1.8 keV to 1 keV as S a increased from 1.0 to 1.4, which also seems 
unphysical, and neither kT e nor r changed their behavior at the interesting values of S a . The new 
model was unable to reproduce the observed counts at energies above 30-40 keV in the island state. 
This, together with the clear change in the behavior of our BB parameters at an "interesting" value 
of S a (or L) seen in Figure 7, leads us to conclude that the BB component parameters listed in 
Tables 3 and 4 are the best fits to our data for the CPL + BB and CompTT + BB models. We 
note, however, that properly measuring low energy components depends on having the correct Nh , 
which we cannot measure with the PCA, and so it is entirely possible that additional components 
exist in the true spectrum. This notion is perhaps supported by the fact that the seed photon 
temperature found in the CompTT + BB model is in good agreement with the BB temperature 
found by BeppoSAX. 
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5.2. Physical Interpretation of the CompTT + BB and CPL + DBB Models 

Our RXTE data on 4U 1820-30 may be described equally well by two different physically- 
motivated models (CompTT + BB and CPL + DBB), both of which have parameters that appear 
to "know" about the saturation of kHz QPO frequency at the accretion rate corresponding to 
S a = 1.3. The physical interpretation of these two models leads to two different pictures of the 
origin of the spectral components in XRBs. Neither picture, however, is completely consistent with 
theoretical and observational constraints. 

The question of the origin of the soft and hard spectral components in XRBs has been addressed 
by many authors (e.g., White, Stella, & Parmar 1988; Guainazzi et al. 1998; Olive et al. 1999, 
Barret et al. 1999, 2000). Early work supposed that the accretion disk intersects the NS surface 
in an optically thick boundary layer which is the source of the blackbody component observed in 
XRBs. A potentially useful diagnostic is the ratio of flux or luminosity in the BB component to 
the flux in the hard component (Lbb / Lh)- Sunyaev & Shakura (1986) showed that, if relativistic 
effects are taken into account, the luminosity of the boundary layer should be at least equal to that 
of the disk, and can be more than twice it if the disk extends all the way to the marginally stable 
orbit. We find Lbb/Lh < 0.5 at all times for the CompTT + BB and CPL + DBB models, and in 
the CompTT + BB model at high accretion rates it is as low as 0.1. Many authors (White, Stella, 
& Parmar 1988; Barret et al. 1999, 2000; Table 1) have found similar results for this and other 
sources. Guainazzi et al. (1998) and Barret et al. (1999) suggest that this implies the weak BB 
component is produced by an optically thick accretion disk, while the strong CompTT component 
originates in an optically thin boundary layer through a process such as the gap accretion model of 
Kluzniak & Wilson (1991). Barret et al. (2000), however, argue that the similarities between NS 
and BH hard X-ray emission require that both be produced by the same mechanism, which rules 
out boundary layer processes. In addition, Church et al. (1998) have found that the evolution of 
spectra during X-ray dips in sources such as 4U 1916-05 are consistent with a compact blackbody 
emitter being occulted while an extended hard X-ray emitter is only partially covered. From timing 
of dip ingresses and egresses they estimate a diameter of ~ 4 x 10 9 cm for the hard X-ray emitting 
region. Thus one possibility is that the BB component is produced in an optically thick boundary 
layer covering a fraction of the neutron star surface, while the CPL or CompTT component comes 
from an extended hot corona. Some mechanism must then reduce the luminosity of the boundary 
layer relative to that of the disk and corona. White, Stella, &: Parmar (1988) suggest that the 
NS may be spinning close to equilibrium with the inner edge of the accretion disk, which would 
lessen the energy dissipated in the boundary layer, but the NS spin periods measured with RXTE 
(typically ~ 300 Hz) are well below this. Another possibility discussed by Barret et al. (2000) is 
that an optically thick boundary layer does not exist. Instead, within the standard optically thick 
accretion disk there exists a hot, optically thin phase of the accretion flow such as the advection- 
dominated solutions of Narayan & Yi (1995). Such a picture has been applied successfully to BH 
systems (Esin et al. 1998). Barret et al. (2000) suggest that the boundary between such a hot 
flow and the NS may well be optically thin and thus neither contribute to the BB component (now 
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coming entirely from the disk) nor quench the hot advective flow and cause it to collapse to a thick 
disk (Yi et al. 1996). The NS surface must still be a source of cool reprocessed photons, however, 
and it is not clear what effect this must have. 

We can interpret our CompTT + BB model fits for 4U 1820-30 in terms of the former of 
these two possibilities, in which the 2.3 keV BB flux comes from an optically think boundary layer 
and the CompTT component from an extended corona. We have noted that in Figure 7 the ratio 
Lbb/Lh decreases as the optical depth of the Comptonizing cloud increases, and both become 
nearly constant at the same accretion rate. In Figure 11 we show Lbb/Lh as a function of r; the 
two are obviously strongly anticorrelated. Thus perhaps the BB emission from the boundary layer 
is simply becoming more obscured as the optical depth of the corona increases, and the boundary 
layer emission is being redistributed to higher energies by the Comptonizing cloud. At the accretion 
rate corresponding to the break, the ~ 2.3 keV boundary layer could be completely obscured, and 
a cooler, ~ 1.7 keV disk component would begin to dominate the BB part of the spectrum. This 
component is emitted from radii outside the inner hot cloud, leading to the upward trend in Rbb 
observed at the same S a . In this scenario a multicolor disk blackbody model would be preferable for 
the higher luminosity spectra, but as noted in Section 3.3 fits of a DBB together with the CompTT 
model were indeterminate. 

Such a scenario ties in with the accretion rate-dependence of the kHz QPO frequencies. The 
origin of kHz QPOs is often explained by means of a "beat-frequency" model. The magnetospheric 
beat-frequency model (Lamb 1991), which was developed to explain horizontal branch oscillations 
(HBOs) in Z sources, has been suggested by Strohmayer et al. (1996) as the source of kHz QPO 
pairs as well. Miller, Lamb, & Psaltis (1998) point out that it is not clear how the X-ray modulation 
at the orbital frequency of the disk is produced in this model, and that it cannot explain both the 
HBOs and the kHz QPOs sometimes seen together in Z sources. Instead, Miller, Lamb, & Psaltis 
(1998) propose the sonic-point beat-frequency model, in which the disk is truncated by radiation 
drag from the NS which causes clumps of gas to lose energy and fall rapidly inward on spiral 
trajectories. This spiral pattern crashes onto the star to create X-ray hot spots which rotate at the 
frequency of the disk inner edge, generating the high frequency QPO. This flow is optically thin, 
which permits the radiation to exert its drag on the disk in the first place. At the same time, some 
gas is channeled by the NS magnetic field onto hot spots which rotate at the spin frequency of the 
star. These hot spots modulate the radiation drag, and thus the brightness of the first hot spots, 
at the frequency of the beat between the disk orbital and stellar spin frequencies, producing the 
second QPO. Miller, Lamb, & Psaltis (1998) point out several predictions of this model as the inner 
edge of the disk, or "sonic point," approaches the marginally stable orbit. First, the frequencies 
of the QPOs should increase together up to a certain accretion rate, then remain constant. This 
is clearly seen in Figure 4. Second, after the marginally stable orbit is reached, the amplitudes 
of the QPOs should decrease. This is because, as M increases further, the optical depth between 
the sonic point and the NS surface will increase, lessening the radiation drag which produces the 
modulation in the first place. Zhang et al. (1998) found that the fractional rms amplitudes of the 
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QPOs steadily decreased as the count rate increased until they disappeared altogether. 

Thus we can envision the following scenario: In the island state the inner edge of the disk lies 
at a radius greater than 6GM/c 2 , within which the sonic-point model is operating. The boundary 
layer is optically thick, producing blackbody radiation from a small equatorial region of the NS 
surface. The NS is surrounded by a hot cloud in the form of an accretion disk corona, since the 
optically thick boundary layer prevents an advection-dominated solution from existing in this case. 
The Wien temperature and radius also indicate that the corona is supported by the disk. The 
electron scattering optical depth of the corona is low, and we may see the boundary layer emission. 
As M increases the inner edge of the disk moves inward and the QPO frequencies increase. The 
optical depth of the corona increases and, as the boundary layer and disk emission increase, the 
corona is cooled and kT e drops. The increased emission from the boundary layer does not make 
it out to the observer, as it goes into cooling the corona. The derived Wien radius also indicates 
that the disk inner radius is shrinking. When the inner edge of the disk reaches the marginally 
stable orbit, the QPO frequencies cease rising. The disk has reached its maximum emitting area 
for cool photons, and so kT e does not fall any more. Further increasing the accretion rate causes 
the optical depth of both the corona and the spiraling gas clumps between the disk edge and the 
boundary layer to rise, reducing the radiation drag that ultimately generates the QPOs. The QPO 
amplitude begins to decrease until, at the point when the boundary layer is completely obscured, 
they disappear altogether. At this point the cooler emission from the disk itself is all that can 
be seen from the outside, resulting in a lower observed UTbb an d higher Rbb, as the boundary 
layer powers the optically thick, nearly thermal-looking Comptonizing cloud. The ratio of BB to 
CompTT flux is quite small at this point. 

The main problem with this picture is that the Comptonizing corona is supported by the 
accretion disk, as indicated by the Wien parameters (Figure 8), and yet its luminosity is always 
more than twice that of the BB component. If the boundary layer is being obscured by the corona, 
transferring energy from the soft to the hard component, then the seed temperature of the CompTT 
model should start to reflect this. Thus if the corona is supported by the disk this picture is not 
consistent with theoretical constraints on the ratio of the flux from the boundary layer and disk 
(Sunyaev & Shakura 1986). 

Our CPL + DBB model for 4U 1820-30 may be interpreted in terms of a hot advective flow 
lying within the standard thin accretion disk if the boundary layer between this hot, optically thin 
flow and the NS may also be optically thin (Barret et al. 2000). The scenario is fairly straight 
forward in this case. In the island state, the inner disk radius is large and the temperature is low. 
Between the disk and the NS the hot advective flow Comptonizes soft photons from the disk and 
NS surface, producing the CPL component. As the disk lies outside the hot corona, it is clearly 
visible. As M increases, the inner disk moves inward and its temperature increases, as shown in 
Figure 9. As the disk area increases, the ratio of DBB to CPL flux increases, as does the cooling 
of the hot flow, causing the cutoff energy to drop. The power law index falls as well, indicating an 
increase in r. The QPOs may be produced by the same sonic point model as before, and so their 
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frequency increases. When the inner disk reaches the last stable orbit, the flux ratio, temperature, 
radius, optical depth, and cutoff energy reach constant values. As M continues to rise, kTsB does 
appear to continue rising slowly as well. 

We may use the DBB parameters at the last stable orbit to estimate the mass of the NS, as 
described by Ebisawa et al. (1994) and Shimura & Takahara (1995). As discussed in Section 3.3, 
Ri n must be corrected by a spectral hardening factor; this effective radius is then set to the radius 
at which the disk temperature becomes a maximum, 5/3 the inner disk radius itself. Then 3Rs = 
(3/5) x 0.6 x Rinf 2 , where Rs is the Schwarzschild radius. The factor of 0.6 accounts for the 
decrease in Ri n \/ cos 9 due to relativistic effects. Using / = 1.7 from Shimura & Takahara (1995) 
and Ri n = 8 km from Figure 9 gives a NS mass of Mjvs = 0.94/a/ cos 9 M Q . The NS mass may also 
be estimated by setting the maximum kHz QPO frequency of ~ 1050 Hz equal to the Keplerian 
frequency at the last stable orbit. This gives Mn$ = 2.1 M for a non-rotating NS. For a NS 
spin frequency of 275 Hz (Zhang et al. 1998) the relativistic corrections are small. Equating these 
two masses requires 9 = 78°, which is nearly problematic due to the lack of dips from this source. 
Should any DBB flux be lost due to scattering in the hot corona, however, the true value of Ri n 
and thus would appear reduced. 

Appealing as this physical model is, it is inconsistent with the observational work of Church 
et al. (1998), who find that X-ray dips in XRBs may be modeled by partial absorption of a large 
corona responsible for the hard X-ray emission. They estimate a diameter of 4 x 10 9 cm for this 
corona, considerably larger than the 16/ 2 = 46 km = 4.6 x 10 6 cm allowed from Figure 9 if the 
hot advective flow is to lie inside the accretion disk. Tomsick, Lapshov, & Kaaret (1998) analyze 
a dip observed from the BH candidate 4U 1630-47 and find that the power law component of the 
spectrum is more highly absorbed, in contrast to the situation found by Church et al. (1998) for 
XRBs. Thus a hot advective Comptonizing flow lying within the accretion disk may be a good 
model for BH systems, but seems to present a geometry problem for XRBs. 

5.3. Lack of Extensive Hard X-ray Emission: NS vs BH Systems 

As noted in Section 2.1, 4U 1820-30 does not belong to the class of XRBs detected at ~ 100 
keV. Our RXTE observations confirm this, as the source is not detected above 50 keV even in the 
hardest state. Extrapolating the island state spectrum from Table 4 into the BATSE energy range 
(20-100 keV) gives a hard X-ray luminosity of ~ 2.5 x 10 36 ergs s _1 , consistent with the upper 
limit of ~ 2.8 x 10 36 ergs s" 1 found by BATSE (Bloser et al. 1996). The predicted flux at 100 
keV in the island state is ~ 2 x 10~ 7 photons cm~ 2 s" 1 keV^ 1 . From Figure 1 is is clear that 
the source spends only a small fraction of its duty cycle in the island state, making hard X-ray 
detection even less likely. Barret, McClintock, h Grindlay (1996; see also Barret et al. 2000) have 
proposed that XRBs can only produce substantial hard X-ray emission when their soft X-ray (1-20 
keV) luminosity is below a critical value of ~ 1.5 x 10 37 ergs s _1 , and that this may be used to 
distinguish NS systems from BH systems. Since the 1-20 keV luminosity of 4U 1820-30 never falls 
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below ~ 2 x 10 ergs s , even in the island state, our lack of detection of the source above 50 
keV is consistent with this picture. 

6. Conclusions 

We have studied the X-ray spectrum of 4U 1820-30 in all three states displayed by atoll 
sources, including the first observation in the island state since the EX OS AT observation described 
by Stella, White, & Priedhorsky (1987). The fact that the 176-day flux variation corresponds 
to normal motion in the CCD indicates that the long period is indeed due to modulation in the 
accretion rate. The spectra are well fit by the CPL + BB, CompTT + BB, and CPL + DBB 
models, but the parameters of the CompTT + BB and CPL + DBB models show more robust 
correlations with each other and with the frequencies of the kHz QPOs when plotted as a function 
of accretion rate (as parameterized by S a ). In the CompTT + BB case, the relative contribution of 
the blackbody component decreases as the optical depth of the Comptonizing cloud increases until 
both reach constant values at the same M at which the QPOs disappear. This can be understood 
in terms of simple geometrical arguments and the sonic-point beat-frequency QPO model of Miller, 
Lamb, & Psaltis (1999) if the blackbody component emanates from an optically thick boundary 
layer and the hard Comptonized component from an extended cloud; however, the relative fluxes 
of the two components violate theoretical constraints. In the CPL + DBB case, the behavior 
of the parameters for the inner accretion disk agree with expectations from the QPOs, but the 
Comptonizing corona must take the form of a hot advective flow within the standard thin disk; 
observations of dipping XRBs suggest the corona is much larger than this. The questions raised by 
both models are central to the issue of hard X-ray production in NS systems. In either case, the 
saturation of QPO frequency at S a = 1.45 provides strong additional evidence for the detection of 
the marginally stable orbit in the accretion disk of 4U 1820-30. 

We would like to thank Dimitrios Psaltis for early discussions of the data. This paper had 
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Fig. 1.— RXTE/ ASM light curve of 4U 1820-30 (2-12 keV) from 1996 October - 1997 October 
(MJD 50400 = 1996 November 12). Each point represents a 1-day average. Pointed PCA/HEXTE 
observation times are marked with vertical lines. The ~ 176 day modulation is clear. 

Fig. 2. — Color-color diagram of 4U 1820-30. Each point represents one of the 90 data segments 
from Zhang et al. (1998), which have typical lengths of 3000 s. The soft and hard colors are defined 
as the ratio of background-subtracted count rates in the bands 3.5-6.4 keV and 2.0-3.5 keV, and 
9.7-16.0 keV and 6.4-9.7 keV, respectively. Filled symbols indicate data that contain kHz QPOs 
from Zhang et al. (1998), open symbols data that do not. The position within the diagram is 
parameterized by the variable S a , the distance along the fitted line. Important values of S a are 
indicated in the Figure for clarity. 

Fig. 3.— Representative PCA (PCU 0) and HEXTE (cluster A) spectra of 4U 1820-30. On 1997 
July 16 the source was in a high state, in the upper banana portion of the color-color diagram. On 
1997 May 2 the source was in the low island state and is detected with HEXTE up to 50 keV. The 
hardening of the spectrum with decreasing luminosity is clear. 

Fig. 4. — Relation between kHz QPO frequency and position within the color-color diagram, 
measured by the parameter S a . Filled and open symbols represent the high and low frequency 
QPOs, respectively. There is a clear saturation in QPO frequency even though S a , and thus 
presumably the accretion rate, continues to rise. 

Fig. 5. — Fitted centroid energy and width (<r) of the gaussian line as a function of total luminosity 
(2-50 keV) for the CPL + BB model. With the exception of the island state observation, the 
centroid energy remains constant at ~ 6.8 keV and the width remains constant at a ~ 0.8 keV, 
indicative of reflection from ionized gas in the disk. The line energy in the island state, ~ 6.55 keV, 
suggests less ionized gas, but the difference in only ~ 2a. A smeared absorption edge at 8.9 keV is 
included in the island state fit. 

Fig. 6. — Examples of 4U 1820-30 spectra fit with the CompTT + BB model. The top panel shows 
the raw count rate spectrum with the fitted model folded through the instruments' responses, the 
middle panel shows the residuals of the fit, and the bottom panel shows the unfolded spectrum 
with the individual model components. The dotted line is the blackbody, the dashed line is the 
CompTT model, and the dash-dot line is the gaussian line. On the left is the island state spectrum 
of 1997 May 2, and on the right is the banana state spectrum of 1997 July 16 (first segment); see 
Table 4 for the fitted parameters. 

Fig. 7. — Variation of spectral parameters for CPL + BB and CompTT + BB models with accretion 
rate, parameterized by S a . Here kTsB = blackbody temperature, Rbb = blackbody radius, a = 
power law photon index, and E c = power law cut-off energy, r = Comptonization optical depth 
and kT e = Comptonizing electron temperature. Error bars are la for one interesting parameter. 

Fig. 8. — The temperature of the seed photon distribution kTy/ from the CompTT + BB model, 
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with the effective Wien radius Rw described in the text. The temperature and radius are anticor- 
related below S a ~ 1.3 and constant above. 

Fig. 9. — Variation of CPL + DBB parameters with S a . Here kT in = temperature of inner disk 
and R in (cos6) 1 / 2 = inner disk radius with inclination 9, while a and E c are as in Figure 7. Error 
bars are la for one interesting parameter. 

Fig. 10. — Comparison of RXTE spectral parameters with previous measurements (see Table 1) 
as a function of total 2-50 keV luminosity. The luminosities of the previous measurements have 
been translated into the 2-50 keV band based on the reported parameters and fluxes. EXOSAT 
= open circle; Ginga = square; ASCA = diamond; Einstein = cross; BeppoSAX = hourglass. The 
luminosities of the RXTE fits have been divided by a factor of 1.15 based on fits to the Crab (see 
text). 

Fig. 11. — The ratio of blackbody luminosity to CompTT luminosity Lbb/Lh as a function of 
r in the CompTT + BB model. The two are strongly anticorrelated; as the optical depth of the 
Comptonizing cloud increases, the relative contribution of the blackbody decreases. 
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Table 1. Previously-Reported Spectral Fits of 4U 1820-30 3 



Date(s) 


N H h 
(10 22 cm" 2 ) 


kT BB c 
(keV) 


a d E c e 


kTj 
(keV) 


r g 
(keV) 


Ei h 
(keV) 


L 










ANS HXX 


(1-28 keV) 








1975 Mar 24-8 

1975 Sep 26-30 

1976 Mar 23-7 


1.3 ±0.3 

< 1.0 
0.93 ± 0.3 




2.12 ±0.08 
1.41 ±0.30 
1.89 ±0.08 






~ 6.7 
~ 6.7 
~ 6.7 


14 
3. 
14 








EXO SAT ME 


(1-30 keV) 









1984 Sep 26 


0.4 ±0.1 


2.0 


1.7 12.0 




6. 


1985 Apr 16 












1985 Sep 22-3 














0.39 


1.61 


3.3 11.9 


6.7 


5. 


1985 Aug 19-20 


0.4 ±0.1 


2.3 


2.5 >30 




2. 




0.52 


2.06 


3.5 13.2 


6.7 


2. 


Ginga LAC (1-20 keV) 


1987 May 1-2 


0.95 ± 0.22 


2.35 ±0.12 


2.21 ±0.14 4.71 ±0.35 


6.7 


5. 


1987 May 4-5 


0.99 ±0.25 


2.43 ± 0.08 


2.16 ±0.16 4.69 ±0.38 


6.7 


6. 


ASCA GIS (0.6-11 keV) 


1993 Oct 6 




0.76 ± 0.02 


3.57 ±0.24 13.5 ±0.6 


6.7 


3. 


Einstein SSS±MPC (0.5-20 keV) 


1978 Apr 7 


0.29 ± 0.002 




1.46 ±0.03 11.4 ±0.2 




5.; 




0.27 


0.89 


3.1 13.3 ±0.7 




5.: 



BeppoSAX LECS±MECS±HPGSPC±PDS (0.3-40 keV) 



1998 Apr 17-18 
1998 Sep 19-20 



0.1 



0.47 



0.55 



4.5 



0.28 ± 0.03 0.46-0.66 



2.83 ±0.08 13.7 ±0.5 



a Only shown are fits using a cut-off power law or Comptonization (CompST in XSPEC) model, plus an option 

b Hydrogen column density using cross sections of Morrison & McCammon (1983) 

c Blackbody temperature 

d Cut-off power law photon index 

e Cut-off power law cut-off energy 

f Comptonizing electron temperature 

g Comptonizing cloud optical depth 
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Table 2. RXTE Observations of 4U 1820-30 





Starting Time 


IVlcan Count Rate 




/'MF r n a 

l IVIHj L 1 


V^lb D 1 V l yz* 1U KcV J 


1996 Oct 26 


88954464 


434.3 


1996 Oct 26 


88977504 


388.5 


1996 Oct 28 


89133104 


481.4 


1996 Oct 28 


89161904 


498.6 


1996 Oct 30 


89294096 


655.2 


1996 Oct 30 


89322896 


624.1 


1996 Oct 30 


89350464 


551.2 


1997 Feb 9 


98100800 


694.6 


1997 Mar 6 


100306640 


785.8 


1997 Mar 7 


100312400 


813.8 


1997 Mar 21 


101575104 


540.6 


1997 Mar 23 


101753664 


623.4 


1997 Apr 10 


103308800 


491.3 


1997 May 2 


105199312 


312.0 


1997 May 10 


105887120 


616.3 


1997 May 28 


107432400 


887.2 


1997 Jun 23 


109688672 


855.2 


1997 Jun 23 


109707200 


851.3 


1997 Jul 16 


111647472 


819.9 


1997 Jul 16 


111670640 


698.5 


1997 Aug 15 


114259928 


752.7 


1997 Sep 10 


116477432 


543.0 


1997 Sep 10 


116488824 


562.5 



a Mission Elapsed Time = seconds since 1994 Jan 1, OhOmOs 
UTC 
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Table 3. Spectral Fits of 4U 1820-30 with the CPL + BB Model 



S a 


Int Time (s) 


kT BB 




Rbb^ 


a 


E c 




Eqw b 




Lbb/Lh 


X 




PCA 


HEXTE 


(keV) 




(km) 




(keV) 


(keV) 


(eV) 








l.Uo 


8481 


2803 


o r ? +0.04 

z - ' -0.04 


U - 88 -0.03 


9 ncr+0.05 
z - uo -0.04 


90 (17+2.91 
zo - u '-2.36 


fi c^Q+0-12 


80 


9 qi 

Z.Ol 


n 9fi 
u.zo 


n 
u. 


191 
1.Z1 


4195 


1362 


9 07+O.O5 
z -°'-0.03 


1 


1 7+0.04 
1 ' -0.05 


9 nn+ - 06 

z - UU -0.06 


11 cc+1.20 


fi 09+o.u 

b - 8z -0.15 


7fi 
1 


9 fi*? 
z.oo 


n 90 
u.zy 


n 
u. 


1 91 
l.Zo 


5908 


1874 


9 41 +0.05 
z - 4i -0.03 


1 


9C+0.04 
zd -0.04 


1 Cfi+0.05 


a /ia+0-63 


« 04+O.I6 
D.84_ 16 


70 


9 7/1 
Z. / 4 


n in 
u.ou 


n 
u. 


i 9^; 

l.ZO 


3519 


1229 


9 4n +0.04 
z - 4u -0.04 


1 


9/1 +0.04 
z4 -0.05 


1 84+0-06 

i -° 4 -0.07 


9 rr+0.73 
8.ZO_ 075 


fi S7+0-18 
b -°'-0.17 


^7 
O 1 


9 8/1 
Z.o4 


n 11 

U.Ol 


n 
u. 


1 98 


9654 


3174 


9 4R+0.04 
z - 48 -0.04 


1 


99+0.03 
zz -0.03 


1 C4+O.O6 
i - b4 -0.06 


c KS+0.44 
O-D8_ 044 


qq+0.18 
u - yo -0.18 


OO 


9 19 
0. 1Z 


n 90 
u.zy 



u. 


1 in 

l.OU 


13763 


4268 


9 4 i+0.04 
z - 4i -0.02 


1 


91 +0.04 
z± -0.03 


1 crc+0.05 
i - oo -0.05 


fi 19+0-31 
b - iz -0.30 


fi 70+O.I8 
b -' y -0.18 


97 


1 99 
O.ZZ 


n 9fi 

U.ZD 



u. 


1 Q9 
l.OZ 


15425 


5006 


9 ci +0.05 
z - di -0.04 


1 


1 7+0.03 
-0.03 


1 40+O.O5 


c 9 c:+0.36 
°- zo -0.37 


ft QO+0.18 
D.88_ 17 


/I K 
40 


Q 9Q 
O.Zo 


n 97 

U.Z ( 


n 
u. 


i i^ 

1.O0 


9938 


2982 


9 c-7+0.05 
z - d '-0.02 


1 


1 q+0.03 

io -0.03 


1 09+0.03 
1 - oz -0.06 


4 - 8b -0.49 


a oq+0.16 


OO 


O.00 


n 9fi 
u.zo 


n 
u. 


1 QS 
l.OO 


6985 


1268 


9 co+0.05 
z -°°-0.06 


1 


m+o.05 

lu -0.03 


1 90+O.O7 
i - z8 -0.11 


4 R7+0.47 
4 - b7 -0.65 


fi 88+ - 15 
b - 88 -0.15 


c;o 
ou 


Q ^8 
O.OO 


n 97 

U.Z ( 


n 
u. 


1 40 


10203 


3130 


9 tr fi +0.05 
Z - OD -0.05 


1 


1 c+0.03 
io -0.03 


1 97+O.O5 


4 qn+ - 31 
4 - yu -0.36 


fi 78+0- 16 
D -'°-0.15 


ou 


1 77 
0. 1 / 


D 94 



u. 


1 /I9 
1.4Z 


3333 


1029 


9 40+O.O7 
z -^°-0.07 


1 


1 Q+0.05 

lo -0.02 


1 04+O.O6 
i - ri4 -0.06 


r c-9+0.37 
°- oz -0.40 


fi oc+0.15 
b - 85 -0.15 


OO 


1 11 
o.ol 


n 99 
u.zz 


n 
u. 


1 4^ 


2546 


855 


9 fi r;+0.04 
Z - DO -0.08 


1 


1 1+0.18 
i± -0.05 


1 i4 +0 - 09 
1 - 14 -0.19 


4 Qfi+0.30 
4.oD_ 1 02 


fi 90+ 016 
D - yu -0.16 


41 
^±o 


1 7q 

O. / a 


n 9^ 



u. 


1 ^1 
1.01 


1840 


604 


9 fi c;+0.04 

z - DO -o.io 


1 


1 c+0.21 
±o -0.07 


1 n4+ - 10 

i - U4 -0.18 


4 1 1 +°- 62 


fi 8^+°- 14 
b - 85 -0.15 


OO 


A 90 

4.zy 


O 9/I 
U.Z4 


n 
u. 


1 KA 
1.04 


6145 


1912 


9 ^4+0.07 
z - O4l -0.07 


1 


i n +0.05 
lu -0.04 


1 ic+0.06 
i - io -0.07 


/I sn+0.32 
4.8U_ 040 


fi 8fi+°- 16 
b - 8b -0.16 


R9 


Q R7 
0.0 1 


O 91 
U.Z1 


n 
u. 


1 

1.0D 


10123 


3135 


9 44+O.O6 
Z.44_ 006 


1 


1 r+0.06 
lo -0.05 


1 17+0.05 
Li7 -0.05 


°- iu -0.23 


R 09+O.14 
b - 82 -0.14 


^0 

ou 


/I 1 
4. 1U 


OIQ 

u. iy 


n 
u. 


1 ^7 
1.0 / 


7844 


2301 


9 c:o+0.07 
z - oo -0.07 


1 


1 1 +0.05 
ii -0.04 


1 19+O.O6 
i - iz -0.06 


4 09+O.29 
4 -° z -0.33 


fi on+0-15 

D.8U_ 17 


^8 
Oo 


/I in 

4. iy 


O 90 

u.zu 


n 
u. 


1 fifl 


5008 


1638 


9 40+O.O8 
z - 4 °-0.07 


1 


nq +0.07 
uy -0.06 


1 19+0.05 
i - iz -0.05 


5 03+ - 24 


fi 8^+°- 15 
D -°°-0.15 


4^1 


4 14 


n 1 7 



u. 


1 fi"? 
l.OO 


566 


184 


9 fil+0-02 


1 


nq+0.08 
uy -0.07 


u - y °-0.01 


4 Q1+0.43 
4 - 4i -0.01 


09+O.I5 
b -° z -0.14 


40 


/I 70 

4. / y 


O 90 

u.zu 


n 
u. 


1 fi^ 
1.D0 


5417 


1668 


9 C9+O.O6 
z - uz -0.09 


1 


1 a+0 12 
1 4-ru. ±^ 

14 -0.05 


U - y8 -0.07 


4 94+O.39 


a nq+0.19 


71 
1 1 


A 1 /I 
4.14 


O 99 
U.ZZ 



u. 


1 R7 
l.D 1 


15967 


5063 


9 4fi+ - 06 
z - 4b -0.05 


1 


9O+0.06 
zu -0.05 


1 n4+ - 05 

i - U4 -0.05 


4 fiX+ - 18 
4 - b8 -0.19 


a 74+O.12 
b -' 4 -0.13 


^1 
01 


/I 71 
4. / 1 


1Q 

u. iy 


n 
u. 


1 7n 
1. 1 u 


9511 


2962 


9 49+O.O7 
z -^ z -0.06 


1 


I 7+0.07 

II -0.06 


1 n o+0.05 


4 09+O.19 
4 - 8z -0.19 


R 7 c;+0.15 
b -' 5 -0.16 


00 


/I 78 
4+0 


17 
U. 1 ( 


n 
u. 


1 71 
i . i o 


9645 


3200 


9 c 4 +0.08 
z -°^-0.07 


1 


1 7+0.04 
11 -0.04 


n 04+O.O5 
U - y4 -0.09 


4 9R+0.25 
4 - zb -0.43 


fi 89+0- 15 
b - 8z -0.16 


74 


4 4^1 


n 91 



u. 


1.75 


3363 


1030 


2-56+Ho 

— U.1U 


1 


90+0.05 
zu -0.05 


0.94+ , 07 

— U. lo 


4.24+0- 7 9 

— U. ( 


6.92+0-is 


77 


4.06 


0.23 


0. 


1.77 


8527 


2708 


9 /(O+0.09 
z - 48 -0.08 


1 


nq+0.08 
uy -0.06 


92+ - 05 

u ' yz -0.05 


4 47+0.20 

-0.23 


6 77+0- 12 

O-l '-0.14 


69 


4.85 


0.16 


0. 


1.80 


5726 


1812 


9 cq+0.09 
z - d °-0.09 


1 


1/1+0.06 
14 -0.05 


n 07+0.06 


4 21+0- 26 

^■ zl -0.40 




90 


4.61 


0.19 


0. 


1.82 


10348 


3097 


9 rrn+0.07 
z - ou -0.07 


1 


9/1+0.06 
Z4 -0.05 


n qi+0.05 
u - yl -0.06 


4 -1Z.+0.22 


6.8210+5 


60 


4.90 


0.21 


0. 


1.85 


6401 


1964 


9 4C+O.O8 
z - 4O -0.08 


1 


9/1+0.08 
z ^-0.06 


U - 88 -0.06 


4 40+0-20 

4 - 4U -0.25 


6 76+0+ 5 

D - ' D -0.18 


68 


5.10 


0.19 


0. 


1.87 c 


466 


161 


9 44+ - 15 

z -^-0.14 


1 


a 7+O.I6 
^'-0.11 


89+ - 09 

u -° y -0.13 


4 41+0-42 
4 - 41 -0.72 


b -sol8:88 


47 


5.54 


0.20 


0. 


1.91 


3327 


1001 


9 c-9+0.06 
z - dz -0.05 


1 


c;i+0.05 
J± -0.04 


1 05+ - 07 

1 - uo -0.07 


4 fi4+ - 33 


6.781°;^ 


43 


5.25 


0.30 


0. 


1.93 


3347 


1023 


9 ci +0.05 
z - oi -0.05 


1 


^4+0.05 
O4 -0.05 


1-0318:8? 


4-68±8:S2 


6 73+0-1 5 

°-' J -0.15 


32 


5.29 


0.29 


0. 


1.96 


3502 


1015 


9 R 9+0.03 
z - uz -0.05 


1 


cn+0.15 
Ou -0.06 


0.86lg;» 


q 00+0.47 

^• 88 -0.70 


fi on+0.00 
b - 8U -o.oo 


51 


5.39 


0.33 


0. 


2.00 


3510 


1037 


2 59+ - 04 
z - oy -0.06 


1 


ci +0.40 
u± -0.04 




4.201™ 


6 89+ - 2 i 

°-° y -0.20 


32 


5.29 


0.37 


0. 



a Blackbody radius 

b Gaussian line equivalent width 

c Total luminosity, xlO 37 ergs s _1 , for a distance of 6.4 kpc, 2-50 keV. This is the raw luminosity derived from 
fit. 

d Smeared edge included at 8.90lg;^ keV 
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Table 4. Spectral Fits of 4U 1820-30 with the CompTT + BB Model 
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a Temperature of seed photons (assumed to follow a Wien law) 
b Comptonization parameter y = AkT e r 2 /m e c 2 
Equivalent Wien radius of seed photons (see text) 

d Total luminosity, xlO 37 ergs s _1 , for a distance of 6.4 kpc, 2-50 keV. This is the raw luminosity derived from 
fit. 
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Table 5. Spectral Fits of 4U 1820-30 with the CPL + DBB Model 
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